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Abstract 

Nb-substituting Fe impurity in KNbC>3 is studied in first-principles supercell 
calculations by the linear muffin-tin orbital method. Possible ways to account 
for the impurity charge compensation are discussed. Calculations are done in 
the local density approximation (LDA) and, for better description of Coulomb 
correlation effects within the localized impurity states, also in the LDA+C7 
scheme. The achievements and problems encountered in both approaches are 
analyzed. It is found that the impurity possess either a low-spin configura- 
tion (with or 1 compensating electron), or a high-spin configuration (with 
2 or 3 compensating electrons), the latter two apparently corresponding to 
practically relevant rechargeable impurity states. 
Keywords: Ferroelectricity , doping, photorefractive effect. 

I. INTRODUCTION 

Impurities in ferroelectric crystals play an important role in many physical applications, 
particularly in what regards optical absorption and the change of dielectric properties under 
illumination, the so-called photorefractive effects. For a review on this subject see, e.g., 
Ref. n. The explanation of the photorefractive effect is essentially based on the assumption 
of the existence of several charge state of an impurity, which can be switched by the drift 
of electrons or holes in the process of illumination. In what regards specifically Fe-doped 
KNb0 3 , one can address Ref. for the description of the models available, and Ref. [||] for 
practical aspects of tuning photorefractive properties. Therefore, much experimental effort 
has been concentrated on the elucidating the electron configuration of impurities in question, 
related lattice relaxation, or other defects coming along with the substitutional impurities. 
Important information of this kind can be extracted from electron paramagnetic resonance 
measurements. Earlier studies for Fe in KNbC>3 [f|] established the presence of Fe 3+ centers 
without local charge compensation and also of other centers, that were interpreted as Fe 3+ 
with local charge compensation, presumably Fe at Nb site with an oxygen vacancy in the 
nearest neighborhood. The subsequent studies 0] included more detailed comparison of the 
spectra with the description based on different structure distortion models. It was suggested 
that the Fe 3+ ion is displaced by ~0.2A towards the neighboring O vacancy. There seems to 
be now general agreement between the experimentalists that Fe preferentially enters the Nb 
site in KNbC>3. The calculations of the energy balance with empirical interaction potential 
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(see Ref. |J], with the potentials from Ref. 0) show beyond doubt that the Fe substitution at 
the K site is energetically unfavorable. In contrary to this, Ref. || claimed that Fe mostly 
substitutes K in KNb0 3 , based on the similarity of angular yield particle-induced X-ray 
emission profiles for K and Fe. This similarity however does not seem highly convincing. 
One should note that the Fe doping in KTaC>3, that is equistructural and equielectronic to 
KNbOs, favors another pattern, with higher probability of entering the K site ||. 

In the situation when the information on the electronic structure of doped systems is 
available from experiment in a rather indirect way, there is a demand for ab initio or other 
reliable calculations which would provide the energy positions of impurity levels in the 
gap and to compare energetically different mechanisms for the charge compensation and 
the relaxation pattern around impurity. To out knowledge, just two attempts have been 
undertaken in this direction up to now. The result of a shell-model calculation by Exner || 
for a Fe at Nb site and an oxygen vacancy as its nearest neighbor is that Fe is displaced by 
0.38A from the vacancy, i.e., contrary to what was supposed based on the analysis of the 
electron spin resonance data ||. Also, the displacement of four "equatorial" O neighbors 
to Fe is 0.19A outwards in the shell- model calculation 0, contrary to the estimated O.lA 
inwards in Ref. ||. The ab initio calculation with the linearized muffin-tin orbital method 
(LMTO) in the atomic sphere approximation (ASA) for one Nb-substituting Fe in the 2x2x2 
supercell of KNbOs, without any oxygen vacancies ||, has shown that the "breathing" 
relaxation of the Oq octahedron is ~0.05A outwards. The shell model calculation for the 
same geometry (also in Ref. 0) as well predicts the outward expansion of 06 of comparable 
magnitude. 

The deficiencies of these previous calculations are that the shell model depends on the 
empirical interaction parameters fitted into it, and normally discards the anharmonic effects, 
and the LMTO-ASA calculation is normally too crude for reliable estimates of lattice relax- 
ation. Moreover, no data on the electronic structure of impurity was published in Ref. 0. 
The aim of our present study is to provide this information, concentrating specifically on 
the treatment of different charge configurations within the localized (i-shell. In particular, 
we compare the treatment within the local density approximation (LDA) with the so-called 
LDA+£7 approach. The calculations are done by a LMTO-ASA method for undistorted 
lattice, without introducing oxygen vacancies. The accurate analysis of lattice relaxation, 
with and without vacancies, needs the use of a full-potential (FP, i.e., with the potential of 
general shape). Such FP-LMTO calculations are in progress and will be reported elsewhere. 



II. CALCULATION METHOD AND SETUP 



The calculations have been done with the LMTO-ASA method 10 



LMTO code of the O. K. Andersen group in Stuttgart The LDA 



in part using the TB- 
-U calculation scheme 

(discussed below) has been implemented in another LMTO-ASA code by V. I. Anisimov. 
In all calculations, we considered the cubic high-temperature phase of KNbOs with lattice 
constant 4.015A; atomic sphere radii were 1.639A for Fe and Nb, 1.964A for K and 1.050A 
for O. Such choice of radii produces the band structure in reasonable agreement with that 
from FP-LMTO calculations ||12[ . According to our experience, the use of additional empty 
spheres in the interstitials of the perovskite structure does not bring any improvement in 
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the LMTO-ASA calculation. In order to account for the lifting of orbital degeneracy at the 
Fe site, the minimal point group symmetry, including only inversion, was assumed in the 
calculations. 

Although KNDO3 is not a typically ionic compound, and, moreover, the covalency of 
the Nb-0 bond is quite important for giving rise to ferroelectric instability, it is useful to 
refer to nominal ionic charges in order to understand why the problem of charge compen- 
sation occurs, when one considers the Fe impurity. Fe is either Fe 2+ or Fe 3+ in oxides; the 
substitution of (nominally) Nb 5+ with Fe 3+ leaves 2 extra electrons in the system, which 
can be removed along with one 2 ~ ion, thus restoring the neutrality. This is one possible 
mechanism of local charge compensation; another, in principle, possible one is that two Fe 3+ 
substitute Nb 5+ and K + . In terms of neutral atoms building the crystal, 30 electrons are 
provided in total in the valence band per unit cell of KNbC>3 (including K3p, which are 
treated as valence band states). If Fe atom with its 8 electrons substitutes Nb, only 5 elec- 
trons can be absorbed by the valence band, leaving Fe nominally in the 3d 3 configuration. 
This is possible, but expectedly highly instable. In order to restore a typical for Fe 3c? 5 or 3d 6 
configuration, one should either provide extra electrons, or remove one oxygen atom, thus 
reducing the capacity of the valence band. There are essentially only two ways to account 
for a charge compensation in the calculation - one either specifies explicitly the configura- 
tion of impurities/vacancies that provides such compensation, or adds extra electrons to the 
system, implying that their donors are in some distant parts of crystal and do not affect the 
local electronic structure. The second way is probably technically easier, because one needs 
only to search for the Fermi energy corresponding to the specified number of extra electrons 
in each iteration. In order to keep the supercell neutral, a compensating positive charge is 
added in the background. 

Our initial calculations did not consider the charge compensation in any way, simply 
treating the neutral supercell with one Nb atom substituted with Fe. As a smallest possible 
supercell for the modeling of isolated impurity, we considered the perovskite cell doubled in 
all three directions, i.e., including 40 atoms in total. The density of states (DOS) resulting 
from this calculation reveals the hybridization of Fe3<i states with those of Nb and O over the 
whole width of the valence band (see Fig.|l]). Moreover, several groups of localized impurity 
states are split off to form flat bands just above the Fermi level. Majority-spin t 2 g states are 
completely filled, whereas all other d-orbitals are roughly half-filled. Magnetic moment at 
Fe is 1.69 /ig and exclusively related to the 3d shell. 

It may be noted that the ti g states remain relatively more localized than e g , as should be 
expected already from fact that the spatial distribution of the former dominates along the 
line directed towards the interstitial between adjacent K atoms. On the contrary, e g states 
directly overlap with the p shells of nearest oxygen atoms. This hybridization mediates the 
interaction between the e g -impurity levels (for both spin directions) at adjacent Fe atoms 
which are situated only two lattice spacings apart. As a result, the e 9 -related subbands are 
smeared so strongly that the band gap almost disappears. We assume that the 8xKNb0 3 - 
supercell is too small for a realistic modeling of isolated Fe impurities. 

The doubling of this supercell with the translation vectors (022), (202), (220) results 
in the \[2 times increased distance between Fe atoms. The broadening of the e g states is 
then sufficiently decreased to form completely split-off discrete levels in the band gap. The 
system remains non-metallic, in contrast to the 40-atom supercell where a small DOS at 
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FIG. 1. Total and Fe3ii-partial DOS as calculated for 40-at. and 80-at. supercells of KNbO-3 
with single substitutional Fe impurity. 
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the Fermi level was present. Therefore, this larger 80-atom supercell seems already to be 
appropriate for the treatment of impurity problem. It should be noted however that the 
charge and magnetic moment at the Fe site are almost identical in the calculation with both 
supercells. 

The number of electrons and magnetic moments within atomic spheres of Fe and its near- 
est neighbors of different types are shown in Table | for several numbers of extra electrons. 
With no extra electrons, the charge configuration of Fe is of course far from the nominal 
3d 3 as mentioned above in reference to a non-charge compensated case. The reason is, the 
electrons of the valence band which have been in perfect crystal localized at the Nb site, 
should be, roughly speaking, added to this number. Since KNb0 3 is not so strongly ionic, 
and because of our choice of atomic sphere radii, the electron numbers in Table | are close 
to those of neutral atoms. Magnetic moment, induced by that of the impurity, primarily 
resides on the first oxygen sphere. 



III. LVA+U TREATMENT 



In a similar supercell calculation for Fe in MgO |JL3|| , it was found that the LDA calculation 
results in a physically wrong stable solution with zero magnetic moment, irrespectively of 
the supercell size. In the present system, Fe impurity has similar nearest neighborhood of 
oxygen octahedron, and the crystal field splits the 3<i-states into lower t 2g and upper e g states 
in a similar way, but the solution is magnetic and stable. Apparently, the difference from 
the MgO:Fe case is in the missing nearest neighbors in the [110] direction in the perovskite 
lattice. In MgO, non-magnetic Mg atoms interact strongly with Fe t 2g states and force the 
occupation numbers in both spin channels to become equal. This mechanism works in the 
LDA, because the potential acting on the t 2g states is not dependent on their occupation, 
and dominates over the intraatomic exchange. 

Although no obvious problem of such kind comes out for for Fe at the Nb site in KNbC>3, 
one should keep in mind that the treatment of localized d states in the LDA may in principle 
be erroneous. A practically feasible way to deal with this problem is the LDA+J7 approach, 



proposed by Anisimov et al. [T3[] and tested by now for many systems, including Fe oxides 



|l5l , |T6H and some perovskites JE7 . The essence of the method is that for the states which are 
ad hoc declared as localized, the total-energy functional as provided by the LDA and casted 
in the form dependent on the total number of localized electrons is corrected depending on 
the actual occupation numbers of different localized orbitals. As a result, the potentials 
acting on localized states become occupation dependent, and the one-electron energies get 
accordingly corrected, thus lifting orbital degeneracies when physically motivated. Effective 
Coulomb interaction U and effective intraatomic exchange J enter as external parameters 
in a self-consistent cycle of LDA+?7 calculation. These parameters can be determined from 
first principles (see, e.g., Ref. once the orbitals which should be treated as localized 

are selected. However, U and J are of essentially intraatomic character and may be largely 
transferred within many systems including the same chemical constituent with a priori 
localized states. 

Solovyev et al. [M thoroughly studied the applicability of LDA+U in comparison with 
LDA to LaMC>3 (M=Ti-Cu) perovskite systems, where a similar problem of treating more 
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FIG. 2. Local 3<i-DOS for charge configurations with 2 and 3 addition electrons per 80-at. 
supercell of KNb03iFe, as calculated in the LDA and LDA+(7. 
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TABLE I. Charges Q and magnetic moments M within atomic spheres of Fe impurity and its 



several neighboring atoms as calculated for different number of added electrons 
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localized t 2g and somehow less localized e g states of a 3d constituent occurs. As apparently 
the most reasonable option (among others tested), the t 2g states have been singled out as 
localized ones, to be treated according the prescription of the LDA+[7 approach, while 
e g have been attributed to common band states well describable within the LDA. This 
segregation may prove useful for the Fe impurity as well. But as a first try, we treated 
all Fe3<i states as localized. Estimates of U for Fe based on somehow different criteria 



vary from ~9 eV |20] to 5.1 eV fL6 |; we used in our calculation [7=6.8 eV as estimated in 



Ref. |n|, and J=0.89 eV. We did the LDA+17 calculations only for the systems with added 
2 and 3 electrons, that correspond to experimentally expected (nominally) Fe 3+ and Fe 2+ 
configurations, correspondingly. The related lines of data in Table | are labeled (+U). 



IV. OVERVIEW OF RESULTS 

As is seen from Table |, the number of electrons on all constituents, at the exception of 
Fe, grows very smoothly with the extra background charge. As the first electron is added, 
nothing happens at the Fe site, and the extra electron goes completely into the valence band, 
resulting only in a slight shift of the Fermi energy but not visibly affecting the DOS. With the 
second extra electron, the Fermi level finally crosses the narrow majority-spin e g subband 
(Fig. ^). This increases the magnetic moment considerably (by less that however, 
because the majority-spin e g states also contribute somehow to the already occupied valence 
band). The distribution of the impurity states in the minority-spin channel is affected only 
quantitatively. The third extra electron begins to populate the minority-spin t 2g subband, 
but most part of it goes into the valence band, where the charge at all sites almost uniformly 
increases. It is noteworthy that the magnetic moment at O neighbors dramatically increases 
with the addition of the second electron, but comes back (and the charge increases, instead) 
with the addition of the third one. Since the states at the top of the valence band are of 



exclusively 02p character (see, e.g., Ref. jTJ]), the drift of the e 9 -subband below the Fermi 
level and into the valence band affects the majority-spin 02p states most directly. The 
crossing of the Fermi level by minority-spin t2 g subband contributes to the minority-spin 
02p subband and reduces the magnetic moment on oxygen atoms, increasing at the same 
time their occupation. At this point, the majority-spin e g subband enters the region of 
noticeable Nb contribution in the valence band and produces the peak of magnetic moment 
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on Nb sites. 

Comparing the results of LDA and LDA+U calculations (Fig. ||]), one should keep in mind 
that the overall effect of the latter is the lowering the energies of occupied states and the 
upward shift of vacant ones. This is exactly what happens in the configuration with 2 extra 
electrons. Since essentially all majority-spin states are already occupied in this configuration 
and all minority-spin states empty, the inclusion of the [/-correction has a negligible effect 
on all integral properties (both charges and moments). However, the exact position of the 
i 2 g-impurity state is changed, that may give rise to the change in the optical absorption. 
The situation in the configuration with 3 extra electrons is completely different. Since here 
the minority-spin t<i g state is partly occupied in the LDA, in the LDA+U treatment the 
potentials acting on the xy, xz and yz components of it become different, thus lifting the 
orbital degeneracy of these states and distorting the local DOS considerably. 

Summarizing, our study of the electronic structure of Fe impurity in the Nb site of KNb03 
with the analysis of different charge compensation has shown, that only two considerably 
different configurations of impurity occur. The first one, with the magnetic moment ~1.7 
/ig, corresponds to the substitutional impurity without compensation - the situation that 
probably is not very common in reality. This configuration, however, survives under addition 
of one extra electron per impurity. Two extra electrons - or, equivalently, Fe impurity 
combined with a distant oxygen vacancy - induces a transition into the high-spin state with 
the magnetic moment ~3.1 /ig, with the minority-spin t 2 g level in the band gap. The charge 
and magnetic moment of this configuration remain intact with the addition of third extra 
electron, but the exact position of the level in the gap (probably, its splitting as well) may 
be affected. The last two configurations with 2 or 3 extra electrons correspond to the most 
practically relevant impurity configurations, referred to as Fe 3+ and Fe 2+ . A more precise 
study of their energetics is possible with the use of a full-potential calculation scheme and 
simultaneous analysis of lattice relaxation around impurity and, optionally, oxygen vacancy. 
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